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ABSTRACT 




Giant resonances in were studied with inelastic 

scattering of electrons with 60 MeV to 105 MeV incident 
energy at a scattering angle of 75 degrees. From 5 MeV 
to 30 MeV excitation energy, five multipole resonances 
were observed of which only the dipole resonances have 
been previously reported. Reduced transition probabilities 
and multipolarity assignments have been made. The five 
resonances were observed at excitation energies of 9.6 
(EO or E2), 11.5 (E2) , 12.3 (El), 15.8 (El), and 23.5 MeV 
(E2) . The splitting of the dipole resonance observed in 
photonuclear reaction experiments with prolate spheroid 
nuclei are confirmed with inelastic electron scattering. 
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I. INTRODUCTION 



Giant resonance investigations using the 120 MeV elec- 
tron linear accelerator o£ the Naval Postgraduate School 
have been conducted since 1973. The previous work has been 
with ^^^Au and ^*^^Pb [Refs. 1, 2, and 3]. The data reduc- 
tion techniques necessary to cope with the large brems- 
strahlung background, called the radiation tail, has been 
greatly improved from the initial attempts. Using computer 
codes and techniques originally developed at Darmstadt and 
reported by Buskirk, et. a7., [Ref. 4] it has been possible 
to determine giant resonances of different multipolarities 
and to make model comparisons with the experimental evidence. 

The previous work at the Naval Postgraduate School 
facility on gold and lead were done at various values of 
momentum transfer, most commonly at a fixed machine energy 
and at varying scattering angles. The investigation even- 
tually covered an excitation range of 5 MeV to 40 MeV. It 
was desired to determine if similar information could be 
achieved by fixing the scattering angle and varying the 
momentum transfer by using different incident energy elec- 
trons. It was decided to study the inelastic electron 
scattering spectra of ^^^Ho in the same energy range for 
three reasons. First, holmium, which is easily obtained in 
a isotopically pure form, has a deformed nucleus with a 
relatively large intrinsic quadrupole moment. Scattering 
from this isotope should illustrate several predicted features 
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unique to deformed nuclear structure. Second, this element 
has been well reported in the literature [Refs. 5, 6, 7, 8, 

9, 10, 11, 12, and 13], particularly the study of giant 
resonances by photoneutron experiments. Since these works 
dealt primarily with the giant dipole resonance, it was 
advantageous to study other multipolarities using the tech- 
nique of inelastic electron scattering. Third, only the 
elastic electron scattering from ^^^Ho has been studied by 

Safrata, et. at. 3 [Ref. 8] and Uhrhane, et. al. 3 [Ref. 12]. 

1 6 S 

Foil samples of ^Ho were prepared and four experiments 
were done at a scattering angle of 75 degrees. Using trans- 
mission geometry, data were collected with incident electrons 
of 60, 75, 90 and 105 MeV energy. The momentum transfer 
squared thus achieved ranged from 0.137 fm’ to 0.420 fm^ . 
Experimental values of the inelastic form factors were de- 
termined from the data and transition multipolarities were 
assigned to five observed giant resonances. Three resonances 
of excitation energy 9.6, 11.5, and 23.5 MeV are reported 
here for the first time. Comparisons between surface oscil- 
lations, Goldhaber-Teller model [Ref. 14] and volume oscil- 
lations, Steinwedel -Jensen model [Refs. 15 and 16] were also 
made . 

There were three objectives of this research. The first 
was to determine if the multipolarities of giant resonances 
could be assigned by measuring inelastic scattering form 
factors as a function of incident electron energy rather 
than the scattering angle as has been done in the past at 
this facility [Refs. 1 and 2]. Second, inelastic electron 
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scattering could be used to investigate multipolarity and 
isospin assignment of the observed giant resonance of the 
^^^Ho nucleus. Third, the splitting of the giant dipole 
resonance is observed in the photoneutron reactions in 
^^^Ho and it was desirable to explore the possibility of 
similar splitting of the quadrupole transitions. 
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II. THEORY 



A. ELECTRON SCATTERING EXPERIMENTATION 

The electromagnetic interaction between the electron and 
the nucleus is well understood and is composed of Coulomb 
interactions between charges and the current and magnetic 
moment interactions of the nucleus with the electromagnetic 
field of the passing electron. All these interactions are 
described by the physical theories of quantum electrodynamics 
and hence are known well enough to make exacting predictions. 
As a consequence, an analysis of electron scattering data 
directly produces information on the target nucleus itself, 
without any details of the imperfectly known strong inter- 
actions interfering in the analysis, such as is the case in 
nuclear structure studies with heavy particle reactions. 

There are two other main possibilities of using electro- 
magnetic interactions for studying nuclear structure, nucle- 
ar absorption and scattering of photons and Coulomb excitation 
by heavy, charged particles. These have disadvantages when 
compared with inelastic electron scattering. Photons cannot 
be utilized to determine ground state ' information since the 
momentum transfer q is determined uniquely by the nuclear 
excitation, o). Electron scattering does not restrict the 
experimentalist in this respect. For electron scattering 
the momentum transfer can be varied through a wide range of 
values since it is given by the expression 
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Cii-1) 




where q is the four vector momentum- energy transfer and 
and k 2 are respectively the four vector incident and scat- 
tered electron momenta. The magnitude of q transferred to 
the nucleus is dependent only on the incident and scattered 
electron energies and the scattering angle by the relation 
[Ref. 16] 

q^ = - 4 sin^(e/2) (11-2) 

for highly relativistic electrons. Since q can be varied 
for constant excitation energy of the nucleus, E^ = E^ - E^, 
momentum analysis of the scattered electrons may reveal 
excitations not previously observed in photonuclear work. 
Coulomb excitation, using the electric field of heavy, 
charged particles, has the same advantages as mentioned 
above for the electron. However, these experiments are 
limited in that magnetic transitions are very difficult to 
excite and the energy of the incident charged particles must 
be limited to less than the Coulomb barrier energy. 

In the first order Born approximation for electron 
scattering the electrons can be considered to be relativis- 
tic Dirac particles scattered from massive point nuclei as 
was done by Mott [Ref. 17]. The incident particle, the 
electron, is assumed to have a spin and a Dirac magnetic 
moment, although the scattering center is assumed to have 
neither spin nor a magnetic moment. The resulting differ- 
ential cross section for scattering is thus given by 
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2 










COS ^ (6/2) 
sin** (6/2) 



(II-3) 



for highly relativistic electrons. For a nucleus of finite 
size a quantity determined by the nuclear structure must be 
included. The measured differential cross section becomes 



(do/dJ^) 



measured 



(II-4) 

(do/dfi)j^j^^^. I/p(?)e^5*^ dt|^. 

nuclear volume 



The factor multiplying the Mott cross section is the "struc- 
ture factor" or "form factor" analogous to electron and 
x-ray diffraction studies, and is in fact called the nuclear 
form factor. It is usually this term which is determined 
in elastic scattering experiments from calculated Mott cross 
sections and measured elastic cross sections. This yields 
the nuclear form factor 



B. INELASTIC ELECTRON SCATTERING 

V 

Analytical expressions for the differential cross section 
for inelastic electron scattering may be derived in a manner 
similar to that for the elastic cross sections using the 
plane wave Born approximations. In this approximation the 
differential cross section (da/dfi)p^^g^ for nuclear excita- 
tion by inelastic electron scattering can be written as a 
sum over the separate cross sections for electric and mag- 
netic multipole transitions. 
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(dcr/df2) E + E 

X X 



CH-6) 



The electric multipole terms are given by Theissen [Ref. 18] 



(d 0 /dQ)g^ = [X(X+1) ' ^B(CX,q)Vg(0) 



- 1 



(II-7) 



+ BCEX,q)Vp.C0)]R 
and those of the magnetic multipoles by a similar expression 



(do/dQ)j^^ = a^ap^q^^k'^B(MX,q)V,j.C0)R' ^ 



(11-8) 



In these equations 

a^ = 47 tX'\x + 1) [(2X + 1)I:]"^ 

k = E /he 
o o 

R = 1 + he (k^/Mc^) (1-COS0) 
with the following definitions of the terms 
X - transition multipolarity 
a - fine structure constant 

E - incident electron energy 
0 - scattering angle 
M - nuclear mass. 

The factor R takes into account the effect of the recoiling 
nucleus on the final state density. The functions Vg(0) 
and V.p(0) account for the longitudinal and traverse compo- 
nents of the four vector nucleus current density resolved 
into components parallel and perpendicular respectively to 
the three vector momentum transfer q. The component parallel 
to q corresponds to an interaction of the electron with the 
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nuclear charge, whereas the perpendicular component is a 
contribution from the nuclear current and magnetization 
densities. For large compared to the electron rest en- 
ergy the functions Vj^C6) and Vrp(0) are defined by Theissen 
[Ref. 18] to be 

V^ce) = ? 2 (l + cos0) / (y-cos0)^ (H-9) 

Vj(9) = %(2y+l-cos0) / (y-cos0) (1-COS0) (II-IO) 

with y = 1 + E ^ / 2E (E - E ) and E the excitation energy. 
For a scattering angle of 75 degrees both Vj^(0) and V,p(0) 

E 

are approximately one. The quantities represent the 

reduced nuclear transition probabilities which will be de- 
scribed in a later section. 

In the plane wave Born approximation (PWBA) the incoming 
and outgoing electrons are described by plane waves. If the 
charge is large such that Za is not much less than unity, 
the electron field is distorted and can no longer be modeled 
by a plane wave. In this case the electron wave functions 
are used which are solutions to the Dirac equation in the 
presence of a field due to the static ground state charge 
distribution of the target nucleus. This model is called 
the distorted wave Born approximation, DWBA, since the ef- 
fect of the nuclear charge on the wave function of the elec- 
tron is to distort both the incident and scattered electron 
wave functions in a manner such that a plane wave solution 
no longer yields satisfactory results. In practice, the 
plane wave solutions of the PWBA are replaced by phase 
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shifted spherical waves which require computer computations 
and tire concept of form factor loses its simple meaning 
[Ref. 19]. 

C. NUCLEAR MODELS 

In this work the models of the nuclear transition charge 
density as introduced by Goldhaber and Teller [Ref. 14] and 
by Steinwedel and Jensen [Ref. 15] were tested against the 
data for comparison and transition strength determinations. 

A brief description of these models and how they are in- 
corporated into the distorted wave Born approximation cal- 
culations for the inelastic form factors is given here. 

The generalized Goldhaber-Teller model as described by 
Uberall [Ref. 16] for any multipolarity assumes the nucleus 
to consist of four interpenetrating fluids, protons with 
either spin up (t) or spin down (4-) and similar states for 
the neutrons. Any two of these states oscillating 180 
degrees out of phase against the other two gives rise to 
three possible collective modes. The original discussion 
by Goldhaber and Teller [Ref. 14] thus corresponds to pro- 
tons, (pf , p4) , oscillating against neutrons, (nf,n4-). 

Other possibilities are (pi ,nt) against (p4-,n4-) and (pt,ni) 
against (p4-,nf). An oscillation of all four fluids in phase 
may produce a fourth, compress ional mode. It requires the 
nuclear matter to be compressible in a monopole vibration. 

In the Goldhaber-Teller model the charge density of the 
ground state PqCv) is assumed to be displaced rigidly such 
that the charge density becomes 
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(II-ll) 



PCr) = PoCr) - d*Vp^(r) 

->• 

where the displacement vector d between the centers o£ the 

neutron and the proton spheres is small. This then de- 

• • 

scribes the dipole motion. Uberall [Ref. 16] shows that 
this may be further generalized to multipole motion if the 
ground state density is assumed to be deformed by a scale 
factor n such that the charge density becomes 



P(r) = p^(r) H- p^^(r) (H-12) 

with the transition charge density 

d p (r) 

= -n r c — ^ ). CII-13) 

The scale factor itself is then expanded in a multipole 
series 

£+k -2 

1 = J “£m (“-I'*) 

ilm 

with = 26^^ (where 6^^ is the Knonecker delta). R is a 
reference radius given to make n dimensionless. A similar 
multipole expansion is possible for the current and magneti- 
zation densities. 

For the Steinwedel and Jensen model a collective motion 
of the neutrons and protons is still assumed but the changes 
are brought about by shifts in the relative densities of the 
two fluids bounded by a rigid nuclear surface of radius R. 

As described by Uberall [Ref. 16] this requires the charge 
density to be modeled by 



P(r) = Pp(t) + P^Cr) = constant (11-15) 
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1/3 

for r < R with R = 1.2 A ' . Again a scale factor is 

-00 

introduced such that 



Pp(r,t) = I pQ + n(r,t) 

" I *^o ■ 



CII-16) 

CII-17) 



As the nuclear surface is assumed to be rigid there can be 
no outflow through the surface hence 



On/9r) = 0 CII-18) 

and again n is expanded in a multipole form. 

These transition densities may then be translated into 
theoretical form factors by using the expression of Ziegler 



[Ref. 19] 



|F(q) 



4tt f X + 1 
Z2 '' X 



(2X+1) 



^ (n-19) 

B(M^,q) 



for the transverse (E) and magnetic (M) components. A simi- 
lar expression can be written for the longitudinal components. 

p 

The coefficients, B(j^X,q), are the reduced nuclear transition 
probabilities mentioned before. Alder, et. al. , [Ref. 20] 
defines the transition probabilities by the relation 



BCX,J.->J^) = |<Jfl |M(X,q)| |J->|" CII-20) 

i 



where the operator M represents the particular transition 
type and and the initial and final angular momentum 
of the nucleus, respectively. Zeigler [Ref. 19] develops 
the representation of these operators using the plane wave 
Born approximation for Coulomb (longitudinal) interactions, 
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M(CX,q) = p^^j^(qr)Y^j^(0,<f)dT, (H-21) 

q 

for the electric part of the transverse interaction 



M(EA,q) = V-L[j^(qr)Y^^(0,^)]dx (II-22) 

and for the magnetic part of the transverse interaction 



M(MA,q) 



. (2A + ir. I 



jn-L[j;^(qr)Yj^j_i(e,<.)]dT CH-23) 



/S 

where is the nuclear current operator as defined by the 
particular model, L is the orbital momentum operator, and 
j^(qr) is the spherical Bessel function of order A. Follow- 
ing the development of Uberall [Ref. 16] the reduced transi- 
tion probabilities become in the plane wave Born approximation. 






_ ? A 



iACqr)PtrCr)dT 



(II-24) 



The model dependent inelastic form factors are calculated 
using the distorted wave Born approximation in the computer 
code GBROW described by Zeigler [Ref. 19] . The reader is 
also referred to the discussion by Uberall [Ref. 16] par- 
ticularly Section 6.4. In this particular code the transi- 
tion probabilities are normalized to unity. To determine 
the strength of the observed transitions it is thus necessary 
to calculate the magnitude of the theoretical form factors 
which would give similar results as the measured form fac- 
tors. The average of the strenghs determined for each 
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experiment is then taken as an estimate of the observed 
transition strengths. The B values thus determined are 
then used to compare the theoretical inelastic form factors 
with the experimental values to determine a best fit to the 
data points. This final B value is then reported as the 
observed strength of the resonance. 

To determine whether or not the observed resonances are 
collective phenomena, comparison of experimental to single- 
particle reduced transition probabilities may be made. The 
ratio of these values would be expected to be significantly 
greater than one for the giant multipole resonances. Further- 
more, a giant resonance should exhaust an appreciable frac- 
tion of the appropriate sum rule. If an observed resonance 
greatly exceeds the sum rule under the assumption of a par- 
ticular transition, it is unlikely that the multipolarity 
assignment is correct. 

The single -particle reduced transition strengths are 
given in Weisskopf units [Ref. 21] 

»CEA)spu = ^ 4-^ ITS R' ' 

and 

BCM^)SPU = — IT 3 " ( 0 . 0111 ) (11-26) 

where R = 1.2 = 6.47 fm for ^^^Ho. 

o 

Another, to some extent model independent, evaluation 
of the transition strength of EX modes may be achieved by 
expressing the strength relative to the energy weighted sum 
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rule (EWSR) . Sum rules for the magnetic transitions are 
more dependent upon the current and magnetic moment distri- 
butions of the nucleus and hence are more model dependent. 
Because of these limitations and since magnetic transition 
strengths are very small at forward angles, they are not 
presented here. The sum rule for isoscalar (AT=0) excita- 
tion modes with L>1 is given by Nathan and Nilsson [Ref. 
22 ], 



S(EX, T=0) = S (E.-E.)B(EX,q) 

f t 1 



Z^e^X(2X+l) ^h^<R^^~^> 
8 IT AM 

P 



(11-27) 



The isovector (AT=0) sum rule for L>1 is related to the 
isoscalar sum rule by 



S(EX,AT=1) = S(EX,AT=0)| . (11-28) 

The corresponding sum rule for an isoscalar monopole (EO) 
excitation is given by Ferrell [Ref. 23] 

1,27 (n- 29 ) 

S(E0) = Z (E^-E.)|M^.|2 = ^ <R^> 

f p 



where is the monopole matrix element. The strength of 

El resonances can be expressed relative to the sum rule 
given by Warburton and Weneser [Ref. 24] 



9e^h^ NZ 

' FFTT — 
P 



(11-30) 



Table I gives values of the energy weighted sum rule (EWSR) 
and the single particle transition strengths for 
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D. GIANT RESONANCE PHENOMENA 



From the study of nuclear level structure by photonucle 
ar processes, it is apparent that there is a collective fea 
ture of the nuclear excitation spectrum called the "giant 
resonance." The term was used initially to describe what 
is now recognized as the giant electric dipole resonance 
(GDR) . In their original paper proposing that these fea- 
tures were indeed resonance structures, Goldhaber and 
Teller by inference also raised the possibilities that 
giant resonances existed with other multipolarities. Ex- 
perimental evidence now indicates that these different 
order multipolarity resonance phenomena do in fact exist. 

Early attempts to explain the giant dipole resonance 
were based mainly on collective models. The initial work 
by Goldhaber and Teller [Ref. 14] mentioned previously as- 
sumed that the protons and neutrons behave as two inter- 
penetrating incompressible fluids. Using a strictly 
classical approach, the two fluids are considered to be 
relatively displaced during dipole oscillations such that 
the proton and neutron fluids no longer overlap near the 
surface. If the restoring force is assumed proportional to 
the surface area, or to R^ , then the frequency of the re- 
sulting harmonic motion is proportional to the square root 

of the force divided by the mass, or w ~ (R /R ) ^~R 
1/3 

Since R~A , then the resulting energy of the harmonic 
oscillator should be E(GDR)~A Goldhaber and Teller 

[Ref. 14] found the relation to be roughly 40A It 
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should be noted here that the original Goldhaber-Teller 
model does not allow for monopole transitions since the two 
fluids are considered incompressible and hence changes in 
relative density cannot occur. The development of the 
generalized Goldhaber-Teller model by Uberall [Ref. 16] 
does make provisions for a monopole transition in the ex- 
pansion of n in equation (11-13). 

The Steinwedel and Jensen model as described by Danos 
[Ref. 25] and in Section II-B assumes a collective motion 
v^ithin a rigid boundary. The energy of the giant dipole 
transition is shown by Danos [Ref. 25] for the spherical 
nucleus to be 



E(GDR) 



2.0 8 fsKNzl 

R M*A2 



(11-31) 



where M* is an effective nucleon mass and K is the energy 
associated \\fith the symmetry energy term of the semi- 
empirical mass formula. Since the nuclear radius is pro- 
portional to then E(GDR)~A for the Steinwedel- Jensen 

model. One selection of the variables as made by Hayward 
[Ref. 26] yields 



E(GDR) 



80 A 



-1/3 



MeV. 



(11-32) 



Since these collective models of the nucleus predict 
that the giant dipole resonance energy varies inversely with 
the radius, a natural consequence is that for deformed 
nuclei having two characteristic dimensions (prolate 
spheroids) the giant resonance should be a superposition of 
two resonances, Okamoto [Ref, 27] and Danos [Ref. 25] have 
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shown that the overall width o£ the dipole resonance is 
strongly correlated to the nuclear deformation as measured 
by the intrinsic quadrupole moment. According to Danos 
[Ref. 25] the connection between the resonance energies E 

3 . 

and E^ and the lengths of the long and short axes a and b 
respectively is given by 



^ = 0.911 ^ + 0.089. (11-33) 

For ^^^Ho the higher energy resonance corresponding to charge 
oscillations along the two short axes of the prolate sphe- 
roid would then comprise two-thirds of the integrated cross 
section. If the width, F, of the total giant dipole reso- 
nance is independent of the energy, Danos [Ref. 25] has 
shown that the ratio of the heights of the separate peaks 
is related by 



) 1 

a'^max _ 

°^b^max ^ 



(11-34) 



Hence, the strengths of the two components of the El 
resonance have the simple ratio 1:2, corresponding to the 
number of degrees of freedom. 

In the earliest work with ^^^Ho by Petree, et. al.^ 
[Ref. 5] a single dipole resonance peak was reported at 
14.5 MeV with a width of 7.5 MeV. Fuller and Hayward 
[Ref. 6] first observed the splitting of the GDR in ^Ho. 
They reported peaks at 12.3 MeV and 16 MeV. Bramblett, et. 
al.i [Ref. 7] found the energies to be 12.10 MeV and 15.75 
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MeV with widths of 2.65 MeV and 4.4 MeV respectively. 

Berg^re^ et. al. , [Ref. 9] reported splitting of the dipole 
resonance at energies of 12.01 MeV and 15.59 MeV. The cor- 
responding widths of the Lorentz shape fits were T = 2.52 
MeV and T = 5.12 MeV. The same splitting as observed by 
Kelley, et. al., [Ref. 10] was 12.32 MeV (T = 2.32 MeV) and 
15.78 MeV (T = 5.04 MeV). With data taken from photoneutron 
interactions with oriented ^^^Ho nuclei, Berman, et. al.^ 
[Ref. 11] reported slightly different energies of 12.28 MeV 
(T = 2.57 MeV) and 15.78 MeV (T = 5.00 MeV). It can be 
shown from their data that the ratio of the observed rela- 
tive transition strengths rather than the apparent inte- 
grated Y cross section of the two resonances is for the 
Saclay group 1.92 [Ref. 9], 2.04 for Kelley, et. at., [Ref. 
10] and for Berman, et . al. ^ [Ref. 11] 1.76. 

Schiff [Ref. 28] developed an analysis of the inelastic 
electron scattering cross section contribution from the 
transition quadrupole moment. Ligensa, et. at. 3 [Ref. 29] 
have extended this work with the Steinwedel- Jensen model to 
include quadrupole multipolarity for deformed nuclei. They 
predict five main quadrupole resonances. The first experi- 
mental evidence of a giant quadrupole (E2) resonance below 
the dipole resonance was found in inelastic electron scat- 
tering experiments at Darmstadt [Ref. 30] . The Sendai group 

20 8 

[Ref. 31] reported a quadrupole resonance in Pb at 22 MoV. 

Those results were confirmed at NPS [Ref. 3] , in experiments 
208 197 

with Pb and Au \'/here E2 resonances were reported at 
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22.5 MeV and 23.0 MeV respectively. The lower E2 resonance 
-1/3 

at E = 63 A ' MeV is generally referred to as the iso- 

-1/3 

scalar (AT=0) and the higher resonance at E = 130 A ' as 
the isovector (AT=1) [Ref. 3] . The existence of the lower 
quadrupole resonance and the observed splitting of the El 
resonance in deformed nuclei makes the evaluation more 
difficult. If one or both of the E2 resonance and the lower 
El resonance are relatively wide and separated by only a 
few MeV then disentanglement of the several resonances may 
be impossible. 

The possibility of a monopole (EO) transition was first 
208 197 

reported in the Pb and Au study mentioned previously 

“1/3 

[Ref. 3] . A mass number dependence of E = 53 A ' was 

X 

found in the data fit analysis. A corresponding transition 
in ^^^Ho has not been previously reported as photonuclear 
absorption cannot excite a monopole transition. 
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Multipolarity 

X 


AT 


Sum Rule 
S(EX,AT) 

eHleV fm^^ 


Single 

Particle 

Strength 

2 ^ 2 X 
e fm 


0 


0 


7.63 X 10“ ^ 




1 


- 




8.10 


1 


1 


5.92 X 10^ 




2 


- 




3.18 X 10^ 


2 


0 


6.16 X 10“ 




2 


1 


9.01 X 10“ 





Units for the monopole transition matrix elements are 
MeV-fm** 



Table I. 



Sum Rules and Single Particle Transition Strengths 



for 
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Ill . DATA ACQUISITION 



A. EXPERIMENTAL PROCEDURES 

It was felt necessary to choose an optimum scattering 
angle which would disentangle the broad resonance structures 
of the giant multipole transitions DWBA calculations of the 
inelastic form factors for El and E2 transitions in ^^^Ho 
were made using the computer code GBROW [Ref. 19] . Both 
volume and surface oscillation models were used for the di- 
pole resonance. Figures 1, 2, and 3 show the results of 
these theoretical form factors as a function of scattering 
angle. It was determined that a forward scattering angle 
of 75 degrees v;ould be optimum for this experiment. At 
backward angles, greater than 90 degrees, transverse mag- 
netic components would contribute significantly to the 
differential inelastic scattering cross section. These 
would complicate this study since it was desired to inves- 
tigate, in particular, the dipole and quadrupole transition 
region of the spectrum. Calculations for the upper dipole 
transition which is due to the oscillation along the shorter 
axes, in ^Ho indicated that this resonance would be suf- 
ficiently large in cross section to be readily detectible 
and not near a diffraction minimum if a scattering angle 
of 75 degrees were used. Also with increasing momentum 
transfer, produced by increasing the incident electron en- 
ergy, the theory predicted a decrease in the relative 
strengths of the dipole resonance as compared to the 
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quadrupole resonance as seen in Figures 1 to 3. The selec- 
tion of 75 degrees as the experimental scattering angle was 

not without some bias. This angle would allow direct corn- 

197 20 8 

parison with the Au and Pb inelastic electron scat- 
tering results reported by Pitthan, et. al.^ [Ref. 3]. 

A foil target of ^^^Ho from Ventron Corporation was 
prepared and roll formed to a thickness of 145 mg/cm^ . The 
target foil was approximately one centimeter wide and three 
centimeters in length. A second foil of 259 mg/cm^ was 
also prepared with similar dimensions. This size permitted 
the same foil sample to be used for all of the experimental 
runs since the focused beam area of the NPS linear accelera 
tor is approximately 0.1 cm^ . The sample was placed in the 
target chamber and positioned for transmission geometry. 
(The electron linear accelerator facility of the Naval Post 
graduate School is described by Warshawsky and Webber [Ref. 
1 ].) 

Four experimental runs were performed at incident elec- 
tron energies of 60, 75, 90, and 105 MeV. All spectra ex- 
cept for 60 MeV investigated excitation energies to 40 MeV. 
The 60 MeV data was taken to 30 MeV excitation energy using 
the 259 mg/cm^ target. The 75 MeV and 90 MeV experiments 
were done using both foils at each energy on separate oc- 
casions. The 105 MeV experiment used only the thinner tar- 
get. It had been determined that using the thinner target 
reduced the background due to the radiative tail as is to 
be expected. In addition, the inelastic spectrum of the 
75 MeV experiment with the 145 mg/cm^ target was observed 
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twice with the same machine parameters to determine if 
there were significant differences in the spectra due to 
machine fluctuations during the course of the 48 hour data 
collection period. These experiments showed reproducible 
results confirming the stability of the experimental con- 
ditions. The count rates were maintained below 20 counts 
per second per channel to avoid losses due to electronic 
dead time of the triple coincidence ladder counting system. 
The total charge delivered by the beam was adjusted for each 
experiment to maintain less than three percent statistical 
uncertainty per channel. Table II lists the experimental 
conditions used for each run. 

B. DATA ANALYSIS 

The line shape fitting procedure developed by Darmstadt 
[Ref. 32] and the radiative tail calculations reported by 
Buskirk, et. at. 3 [Ref. 4] have been incorporated into a 
computer code called NAW [Ref. 33] for an IBM 360/67. This 
code is described and presented in Ferlic and Waddell [Ref. 
2]. The entire inelastic spectrum including the resonances, 
radiation tail, and background were simultaneously fit as- 
suming Breit-Wigner shapes for each resonance. Initial 

resonance energies and widths were chosen based on previous - 

-1/3 

ly reported photoneutron data and from the general A ‘ 
behavior of giant multipole resonances [Ref. 3] . Resonance 
positions and widths were varied as input parameters v\fhile 
the computer code adjusted peak strengths and background 
parameters to obtain a best fit to the experimental data. 
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The code was subsequently modified to provide an option 
for maintaining a specifix: ratio of the two dipole reso- 
nance components. This was done to preserve the known 
giant dipole relative strength ratio as indicated in the 
theory and found by photonuclear experiments. (See Section 
II-D.) 

The criteria used to determine a reasonable fit were 
as follows: 

(1) The data and calculated spectrum should coincide 
visually on the superimposed spectrum plots. 

(2) The degree of freedom should be less than 

one. The data is not strictly statistical because the de- 
tector momentum interval is larger than the momentum in- 
crement of the spectrometer field and hence correlations 
exist between energy bins. 

(3) All observed resonances and widths should consis- 
tently fit each of the four spectra. 

The 60 MeV spectrum is not considered as reliable as 
the other spectra even if these criteria were not apparently 
violated due to the use of the thicker target and the low 
incident electron energy. Here slightly different back- 
ground parameters from those used at other energies were 
necessary to fit the data because of the differing experi- 
mental conditions. The higher background due to the radia- 
tive tail made extraction of the experimental inelastic 
form factors from the data difficult and somewhat arbitrary. 

Figures 4 through 6 represent the experimental inelastic 
spectra corrected for spectrometer dispersion effects with 
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the fitted total background and the individual resonances 
superimposed. Figures 7 through 9 present corresponding 
spectra with the total background subtracted revealing the 
resonance structures and their features. 

The areas under the inelastic resonances were calcu- 
lated from the Breit-Wigner shapes and the ratio of each 
inelastic resonance area to the area of the elastic peak 
were presented in the NAW output with a statistical uncer- 
tainty. The elastic form factors for the experimental 
energies were determined using the phase shifted spherical 
wave calculations of Fisher and Rawitsher [Ref. 34] . The 
inelastic form factor could then be determined by the 
relation 

|Fil" = CAi/Aei)|Fell'- (IH-1) 



For these calculations the ground state charge distribution 
was assumed to be of the Fermi form 



PoCr) ■= 



1 + e 



(r-c/4 . 4t) 



-1 



CIII-2) 



with c = 6.18 fm and t = 2.51 fm from Uhrhane, et. at. 3 
[Ref. 12] . Table III presents the inelastic form factors 
obtained for the 75, 90, and 105 MeV experiments. 

Table IV presents the transition probabilities that 
were obtained with the best fit resonance positions and 
widths. The values are reported for both conditions utilized 
within the computer code NAW - free fitting of all observed 
resonance strengths and free fitting of all strengths ex- 
cept that corresponding to the lower component of the dipole 
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transition. The strength of the latter was maintained at a 
fixed ratio to the upper component. The value of this ratio 
was determined from the (y,n) experiments, using inelastic 
form factor calculations. For the El resonance the experi- 
mental strengths were determined for the surface oscillation 
(Goldhaber-Teller) model and the volume oscillation (Stein- 
wedel- Jensen) model. Comparisons of these results with the 
appropriate energy weighted sum rule and single particle 
transition strengths are also made in Table IV. 

C. ERROR ANALYSIS 

Certain corrections must be made to the observed experi- 
mental electron scattering cross section. The observed area 
of both the elastic and inelastic lines must be increased 
since some of the electrons scattered at the corresponding 
energy suffer energy losses due to radiation during scat- 
tering, radiation and electronic collisions before and after 
scattering, and ionization straggling effects and hence are 
not counted at the proper energy [Ref. 16] . For the NAW 
code the radiation tail is calculated using the method de- 
scribed by Ginsberg and Pratt [Ref. 35] which was modified 
to use a phase shift calculation of the elastic scattering 
cross section to replace the original Born approximation 
cross sections [Ref. 4] . An analytic function is then used 
for the fit of the background due to the radiation tail and 
experimental processes. For this work an expression of the 
form 
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(III-3) 



?! + P 2 A + Pj RDT)exp ^ (E^-Ef) 

o 

was used. is the incident electron energy, E£ the energy 

o£ the scattered electron. A = E - E^ where E^ is an ex- 

o to to 

pansion energy, and RDT is the calculated radiation tail. 

The P^'s are free parameters and are determined by the 
fitting routine. 

Warshawsky and Webber [Ref. 1] describe the fractional 
error associated with the determination of the areas as a 
result of fractional errors in the Poisson counting statis- 
tics, variations in beam intensity during a run, and the 
error associated with the determination of the areas. Ferlic 
and Waddell [Ref. 2] investigated the energy and width re- 
lationships to the error in the determination of the area 
of Breit-Wigner shapes. They found contributions on the 
order of 20 percent in addition to the NAW statistical 
error. The total uncertainty would then be some function 
of these combined effects. 

In addition a correction to the line shape due to a 
change in the inelastic form factor itself over the width 
of the resonance must be considered [Ref. 32] . Calculations 
o£ the theoretical inelastic form factors indicated a change 
in I I ^ of several percent per MeV. For wide resonances 
the energy shift in the apparent resonance peak may be 
significant . 

Repetition of the 75 MeV experiment during the same run 
with the machine parameter unchanged showed that current 
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control techniques maintain experimental conditions over 
the duration of the experiment. Hence, beam variations do 
not seem to be as large now as reported by Warshawsky and 
Webber [Ref. 1] . To provide an estimate of the true frac- 
tional error associated with each resonance would be most 
difficult. For this work an uncertainty of twice the sta- 
tistical error computed by the NAW code is assumed for the 
data. This uncertainty is reported in Table III as a per- 
centage error of the inelastic form factors obtained exper- 
imentally. The true uncertainty is at most this large and 
quite probably less. 
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Experimental 

Incident 

Energy 

(MeV) 


Target 

Thickness 

(mg/cm^) 


Measured 

Incident 

Energy 

(MeV) 


Elastic 

Peak 

Resolution 

(MeV) 


(fm^) 


[Ref. 34] 


60 


259 


60.77 


.40 


.137 


2.6x10" ^ 


75 


145 


76.70 


.42 


.214 


9.41x10"^ 




259 


75.50 


.40 


.214 


9.41x10"^ 


90 


145 


90.85 


.51 


.309 


5.49x10"^ 




259 


90.11 


.53 


.309 


3.49x10"^ 


105 


145 


104.94 


.58 


.420 


1.99x10"^ 



Scattering angle was 75 degrees for all experiments. 
Transmission geometry was used for target position. 



Table II. Experimental Conditions. 



Excitation 

Energy 

(MeV) 


Incident Electron Energy 


7 5 MeV 


1 Error 


90 MeV 


^ Error 


105 MeV 


% Error 


9.6 


1. ISxlO""* 


34 


8.32x10"® 


30 


7.30x10"® 


20 


11.5 


5.92x10"'* 


16 


3.59x10"'* 


12 


2.49x10""* 


12 


12.3 


3. 05X10""* 


38 


9.56x10"® 


10 


6.12x10 ^ 


30 


15.8 


6.62x10""* 


38 


2.45x10""* 


10 


1.53x10""* 


30 


23.5 


4.60x10"“* 


50 


1. 78x10""* 


22 


2.17x10""* 


24 



Fixed Ratio of El Components 



Table III. Inelastic Form Factors. 
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Reduced Transition Probabilities, Fixed Ratio o£ El Components 
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Figure 1. Inelastic form factor for El (upper component) 
-transition predicted by Steinwedel -Jensen . 

Model for 75, 90, and 105 MoV incident electron 
- - - energy, t - i - - 
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Figure 3 



Inelastic form factor for E2 transition predicted 
by, Goldhaber-Teller model for 60, 75, 90, and 
•105 MeV incident electron energy. , j 



.L.,_ 
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Figure 4. 75 MeV inelastic scattering spectrum with back- 

ground . 



43 



Counts per yC 



Typical Uncertainty Shown 



A 




Figure 5. 90 MeV inelastic scattering spectrum with background. 
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Figure 6. 



105 MeV inelastic scattering spectrum with back- 
ground . 
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Figure 7. 75 MeV inelastic spectrum without background. 
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Arbitrary Unit = 1.59x10 '* £m^ / MeV St 




Figure 8. 90 MeV inelastic scattering spectrum without 

background. 



47 



d^a / dEdfi , Arbitrary Units 



Typical Uncertainty Shown 
Arbitrary Unit = 7.31x10 ^ £m^ / MeV St 




Figure 9. 105 MeV inelastic scattering spectrum without 

background . 
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IV. DISCUSSION 



A. COLLECTIVE RESULTS 

The data analysis indicates the presence of five 
resonances at excitation energies of 9.6, 11.5, 12.3, 15.8, 
and 23.5 MeV in the inelastic electron spectrum of 
In the original data which extends to 40 MeV excitation 
energy, there exists the possibility of a sixth resonance 
at 33 MeV. However, it is not certain at this time if the 
peak is a true resonance or is generated by the experimental 
arrangement of the target chamber itself due to elastically 
scattered electrons passing through the metal of the spec- 
trometer collimater. There are two reasons for accepting 
the resonance as real. First, in the energy loss experi- 
ments of Buskirk, et. at., [Ref. 36], it was found that if 
the energy loss in lead is 6 MeV, the width of the distri- 
bution is 10 MeV. If it is now expected that the loss is 
over 30 MeV, the width would be over 30 MeV as compared to 
the 6 MeV width observed in lead. Since the observed struc- 
ture has a preliminary fitted Breit-Wigner shape approximate- 
ly 7 MeV wide, an energy loss mechanism for elastically 
scattered electrons does not seem to be the source. Second, 
an inelastic electron scattering experiment was performed 
with the spectrometer entrance collimator plugged and the 
resulting spectrum had no structure resembling the observed 
33 MeV resonance. The exclusion of a 33 MeV resonance 
would not make significant changes in the transition 
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strengths reported here since the other resonances are all 
lower in excitation energy. 

The conclusion that the other levels are true giant res- 
onances is supported by the comparison of the transition 
strengths with the appropriate sum rule and the single par- 
ticle transition strengths. In Table IV it can be seen 
that an appreciable percentage of the sum rule associated 
with each transition type is accounted for by the observed 
energy weighted B value. In addition, the B values are 
several Weisskopf units indicating that a collective response 
is being observed. 

In Figures 8 and 9 it can be seen that the apparent 
heights of the resonances at 11.5 MeV and 23.5 MeV increase 
relative to the apparent height of the peak at 15.8 MeV as 
the incident electron energy increases. Recall that the 
scattering angle of 75 degrees was selected to observe the 
possible enhancement of the quadrupole resonance over the 
dipole for increasing momentum transfer. Since the 12.3 
MeV and 15.8 MeV resonances are known to be dipole reso- 
nances from the photoneutron work with the first as- 

sessment of the 11.5 MeV and 23.5 MeV resonances is that 
they are quadrupole resonances. 

For each experiment the overall fit was improved by 
the modification to the computer code NAW as described in 
Section III-B. Resolution of the two resonances at 11.5 
MeV and 12.3 MeV was not sufficient due to the large width 
of both peaks. In the initial free fitting attempts it 
was found that the 11.5 MeV resonance was enhanced and 
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that only a relatively small resonance at 12.3 MeV was re- 
quired to meet the criteria of a good fit. This caused the 
calculated B values for the 12.3 MeV resonance to be much 
smaller than the B values of the photoneutron work (Table 
IV) . Also the assumed lower El resonance component had B 
values in Weisskopf units of approximately unity. This 
would indicate a single particle scattering effect. Fit- 
ting with a constraint on the ratio of the strengths of 
the lower and upper dipole resonances gave better agreement 
between measured and calculated form factors. 

B. RESONANCES AND CHARACTERISTICS 

Each of the observed resonances is discussed in this 
section * 

1. 9.6 MeV 

- 1/3 

An EO isoscalar resonance at 53 A was reported 

197 208 

previously at this facility [Ref. 3] in Au and Pb . 

This would correspond to an excitation energy of 9.6 MeV in 
165 

Ho. The existence of a resonance in this region is 
clearly seen in the 105 MeV experimental data in Figure 9. 

The calculated magnitude of the matrix element of an as- 
sumed monopole resonance would exhaust 28 percent of the 
monopole sum rule. 

2. 11.5 MeV 

The resonance at 63 A 11.5 MeV for ^^^Ho, is 

now generally believed to be the isoscalar quadrupole giant 
resonance [Ref. 3] . As noted in Table IV the energy weighted 
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observed strength exhausts 73 percent o£ the quadrupole sum 
rule for the natural width of the resonance assumed. How- 
ever, fitting attempts were made with the width of this E2 
between 2.9 MeV and 4.5 MeV and the obtained in each was 
approximately the same, but with rather different percent- 
ages of the sum rule. The strength of the resonances at 
9.6 MeV and 15.8 MeV were not significantly affected in 
these different iterations. A relatively wide resonance 
would, of course, indicate unresolved quadrupole splitting. 
The uncertainty in the observed width precludes answering 
that question based on this data. 

3. 12.3 MeV 

This is known to be the lower excitation energy of 
the giant dipole resonance for the deformed ^^^Ho nucleus 
[Refs. 9, 10, and 11] . It was found that the excitation 
energy and natural width as reported by Berman, et. al. ^ 

[Ref. 11] in the photonuclear work for the resonance at 
12.3 MeV fit remarkably well to the observed experimental 
data from inelastic electron scattering. This El resonance 
was fitted by a fixed ratio of its inelastic form factor 
to that of the resonance at 15.8 MeV, hence its reduced 
transition strength is a function of the strength observed 
for the upper, better visible, component of the dipole 
resonance. These strengths are discussed in the next section 
An assignment of an El isovector multipolarity based on the 
results of the photoneutron data and of this work is con- 
sidered well supported. 
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4. • 15.8 MeV 



This is the predicted partner in the dipole reso- 
nance splitting of the deformed nucleus ^^^Ho [Refs. 9, 10 
and 11]. The energy and width reported by Berman, et. al.^ 
[Ref. 11] were used in the initial fit and were found to 
not vary significantly in subsequent trials. The experi- 
mental inelastic form factors shown in Figure 10 can be 
described reasonably well with either model. The close 
agreement for the upper component of the dipole resonance 
is considered valid evidence that the observed splitting 
is a real effect in the prolate spheroid, Comparison 

of the energy weighted sums provides some model selectivity. 
The Goldhaber-Teller model exhausts 107 percent of the El 
resonance sum rule with the 12.3 MeV and 15.8 MeV - combined. 
For the Steinwedel -Jensen model the two comprise only 67 
percent of the sum rule. The strengths as reported by the 
Saclay group [Ref. 9] comprise 120 percent of the sum rule, 
while the values reported by Berman, et. aZ., [Ref. 11], 
exhausts 106 percent of the sum. Thus the Goldhaber-Teller 
model is favored by this data in comparison with the results 
of the photoneutron work. The assignment of an El isovector 
to this resonance and the 12.3 MeV resonance is consistent 
with previous work [Refs. 3, 9, 10, and 11]. 

5. 23.5 MeV 

This resonance has not been reported previously. It 
appears strongly in the 105 MeV experimental data. Figure 9. 
Berman, et. al.^ [Ref. 11] noted some structure in the photo- 
neutron work of ^^^Ho in this excitation region. The 
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calculated B value for this work exhausts 49 percent of the 

energy weighted sum rule for an E2 isovector resonance. In 
197 208 

the Au and Pb work [Ref. 3] the isovector quadrupole 

-1/3 

was reported at 133 A ' . This would be approximately 24 

MeV for In addition an E2 assignment would be con- 

sistent with the structure reported by Berman, et. at., 

[Ref. 11]. However, the data because of its relatively 
large width of 7 MeV, indicate a possible splitting of the 
isovector E2 resonance in several levels, which makes a 
fit by one Breit-Wigner resonance quite inadequate. 
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V. CONCLUSIONS 



Giant multipole resonances studied by inelastic electron 
scattering in ^^^Ho indicate five resonances in the excita- 
tion range of 5 MeV to 30 MeV. These resonances occur at 
9.6, 11.5, 12.3, 15.8, and 23.5 MeV. The 12.3 MeV and 15.8 
MeV resonances have been reported in the photoneutron work 
and are known to be electric dipole. In spherical nuclei 
these dipole resonances are threefold degenerate in energy, 
whereas for prolate spheroids such as ^Ho the observed 
single resonance is split into two separately resolved res- 
onances. There exists a correlation between the transition 
strengths of the two resonances which are found to be in 
the ratio of 1:2 for the lower energy resonance in respect 
to the upper component. 

The experiments fit El form factors for either the 
Goldhaber-Teller or Steinwedel -Jensen models. However, only 
the strength calculated from the Goldhaber-Teller model 
agrees with the sum rule and (Y>n) experiments. 

In this series of experiments a significant difference 
in technique compared to most previous experiments consisted 
in the use of one scattering angle and various energies. 

The advantages and disadvantages of this procedure can be 
summaried as follows: 

1. No magnetic transitions appeared in the excitation 
energy range of 5 to 15 MeV, which includes the EO, E2 and 
El transitions. 
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2. The resulting variation o£ the momentum transfer 
q is small, so no diffraction minima are included. If 
nothing were known about these resonances, these experi- 
ments alone would not decide between El and E2 assignments 
as conclusively as the other experiments. 

The giant quadrupole which appears as a general feature in 
heavy nuclei [Ref. 3] is observed at 11.5 MeV for 
One objective was to choose a range of momentum transfer 
such that the E2 resonances would be enhanced relative to 
the El resonance. This enhancement occurred in the 105 MeV 
run, in which the E2 is about twice as large as the lower 
El. The width of the E2 was not decisively determined be- 
cause of the presence of the 9.6 and 12.3 MeV levels. Width 
from 2.9 to 4.0 MeV fit reasonably well, with 4.0 MeV 
favored. Thus the expected splitting or broadening of the 
E2 resonance in the deformed nucleus is not definite. The 
surface oscillation model is strongly favored. The addi- 
tional E2 isovector resonance at 23.5 MeV has not been pre- 
viously reported in ^^^Ho although the general feature was 

reported by Pitthan, et. al. ^ [Ref. 3] in gold and lead at 
-1/3 

133 A ' . An assumed monopole resonance at 9.6 MeV was 

found in the experimental data which corresponds closely 
with the 53 A feature. 

Using the information for the inelastic form factors, 
the corresponding reduced transition probabilities and com- 
parisons with the energy VNreighted sum rules it was possible 
to determine multipolarity of the resonances. Good agreement 
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was found between this experiment and the previous gold and 

- 1/3 

lead experiments in terms of the general A ' relations 
for the various resonances. The excitation energy, width, 
and multipolarity assignments made in this work are pre- 
sented in summary in Table V along with the results reported 
earlier for ^^^Au and 
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Table V. Giant Resonances 



Ex 

A'^^hleV 


Ex 

(HeV) 


^nat 

(MeV) 


Multipolarity 

Assumed 




16 


^Ho 




53 


9.6 


2.0±5. 


EO or E2 isoscalar 


63 


11.5 


4.0±1. 


E2 isoscalar 


67 


12.3^ 


2.5^ 


El isovector 


86 


■15.8^ 


5.0^ 


El isovector 


133 


23.5 


7.0±1 


E2 isovector 




19 


^Au 




53 


9 . 2 


2.2±0.5 


EO isoscalar 


63 


10.8 


2.9±0.2 


E2 isoscalar 


81 


14.0 


4. 5±0.2 


El isovector 


133 


23.0 


7 ± 1 


E2 isovector 




20 


^Pb 




53 


8.9 


1.8±0. 5 


EO isoscalar 


63 


10.5 


2.8±0.3 


E2 isoscalar 


81 


13.6 


3.9±0. 1 


El isovector 


133 


22.5 


5 ± 1 


E2 isovector 



^[Ref. 11] 
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